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Abstract. We have recently examined a large number of points in the parameter 
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these points satisfied existing experimental and theoretical constraints. This analysis 
provides insight into general features of the MSSM without reference to a particular 
SUSY breaking scenario or any other assumptions at the GUT scale. This study opens 
up new possibilities for SUSY phenomenology at colliders as well as in both direct and 
indirect detection searches for dark matter. 
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1. Introduction 

Supersymmetry (SUSY) represents an appealing possibility for Beyond the Standard 
Model Physics; its discovery would help provide answers to many of the preeminent 
questions in particle physics, astrophysics, and cosmology. However, as no sparticles 
have been observed, it is clear that if SUSY exists, it must be broken. The mechanism 
which could break SUSY is a question of great importance, and there is an ever growing 
list of possible scenarios, including mSUGRApJ, GMSB[2], AMSB[3], and gaugino 
mediated supersymmetry breaking[l]. In each of these scenarios, the SUSY spectrum 
is described by a handful of parameters, generally defined at the SUSY breaking 
scale; straightforward RGE running of sparticle masses and coupling constants yields 
predictions for the mass spectra and decay patterns of the various sparticles at energy 
scales relevant for colliders or cosmology. However, these SUSY breaking scenarios are 
restrictive and predict specific phenomenologies for colliders and cosmology that do not 
represent the full range of possible SUSY signatures. 

It is clearly desirable to study the MSSM more broadly without making simplifying 
assumptions at the high scale that may turn out to be unwarranted. However, the MSSM 
requires 105 parameters to describe SUSY breaking, in addition to the parameters of the 
SM|5]. Obviously this is far too many parameters to study directly, so some simplifying 
assumptions must be made. Here we will restrict ourselves to the CP-conserving 
MSSM {i.e., no new phases) with minimal flavor violation (MFV)[6]. Additionally, 
we require that the first two generations of sfermions be degenerate as motivated by 
constraints from flavor physics. We are then left with 19 independent, real, weak-scale, 
SUSY Lagrangian parameters, namely the gaugino masses Mi 2,3, the Higgsino mixing 
parameter /i, the ratio of the Higgs vevs tan/3, the mass of the pseudoscalar Higgs boson 
m^, and the 10 squared masses of the sfermions (m^i 3, m^ji 3,mj^ 35'"^[i 3, and mgi 3). We 
include independent A-terms only for the third generation (Aj,, At, and Ar) due to the 
small Yukawa couplings for the first two generations. This set of 19 parameters has 
been called the phenomenological MSSM (pMSSM)[7]. 

To study the pMSSM, we performed a scan over this 19-dimensional parameter 
space assuming flat priors for the specified ranges [8]: 

100 GeV <mj<l TeV , 
50GeV < |Mi,2,/i| < ITeV, 
100 GeV < M3 < 1 TeV , 

IA*,r|<lTeV, (1) 

1 < tan/5 < 50, 

43.5 GeV < < ITeV. 

The value of 43.5 GeV in the last constraint was chosen to avoid the possible on- 
shell decay Z hA. We randomly generated 10'' points in this parameter space and 
subjected them to a number of existing theoretical and experimental constraints. We 
also performed a scan with log priors and slightly different mass ranges (that we will 
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not employ here) in order to gauge the influence of priors on our results; we found that 
these results are substantially similar to those obtained in our flat prior scan[8]. Using 
these parameters we generate a SUSY spectrum utilizing SuSpect2.34[7]. By convention, 
the parameters are specified at the scale given by the geometric mean of the two stop 
masses. The input values for the SM parameters used in our analysis are given in [8]. 

We then apply a series of constraints obtaining a set of models that satisfy all 
existing theoretical and experimental data. (This is the so-called "fiat prior" set obtained 
in Ref. [8] .) In the analysis below we will discuss the applied constraints then will examine 
the properties of the LSP at the parameter points which remain viable. In particular, we 
will examine the gaugino and Higgsino content of the LSP (which is always the lightest 
neutralino). We will also discuss the nature of the nLSP and the difference between its 
mass and the LSP mass; this is important, for example, in coannihilation processes. We 
will then examine the signatures in direct and indirect WIMP detection experiments 
obtained for these points in parameter space. 

2. Theoretical and Experimental Constraints 

We now discuss the theoretical and experimental constraints which we applied to the 
generated parameter space points (which we shall hereafter refer to as "models" for 
convenience) . We will present each of these briefly in turn; for more details, one should 
consult [8]. 

2.1. Theoretical Constraints 

We demand that the sparticle spectrum not have tachyons or color or charge breaking 
(CCB) minima in the scalar potential. We also require that the Higgs potential be 
bounded from below and that electroweak symmetry breaking be consistent. We assume 
that the LSP, which will be absolutely stable, be a conventional thermal relic so that 
the LSP can be identified as the lightest neutralino. If it is a significant component of 
the dark matter, the LSP must be uncolored and uncharged, thus the LSP can only 
be a sneutrino or a neutralino. The possibility that the LSP is a sneutrino can be 
easily eliminated in the pMSSM by combining several of the experimental constraints, 
particularly those involving direct detection of sneutrino WIMPs and the invisible width 
of the Z, as discussed below. 

2.2. Low Energy Constraints 

The code micrOMEGAs2.20[9] was used to evaluate the following observables for each 
point in the parameter space: Ap, the decay rates for 6 — 57 and Bg and 
the g — 2 oi the muon. In addition, we evaluate the branching fraction for B tu 
following[10j and [11]. The ranges that we allow for these observables are listed in 
Table [H The large range for the SUSY contribution to g — 2 (~ 6a) is due to the 
evolving discrepancy between theory and experiment [12j. 
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Table 1. Ranges allowed for various low energy observables in our analysis. 

Wp implpmpntpd constraints from mpson-antimpson mixing [T3] by assuming 
MFV[6], imposing first and spcond gpnpration mass dpgpnpracy, and dpmanding that 
tlie ratio of first /spcond and third gpnpration squark soft brpaking massps (of a givpn 
flavor and hphcity) diffpr from unity by no morp than a factor of 5. Wp also impospd 
analogous rpstrictions in thp slppton spctor. 

2.3. LEP Constraints 

Wp now consider thp constraints that arise from LEP data. Dup to running LEP at thp 
Z polp, it is vpry unlikply that thprp can bp chargpd sparticlps with massps bplow Mz/2. 
Thp samp constraint is applipd to thp lightpst nputral Higgs boson. Data from LEPII[22] 
suggpsts that thprp arp no npw stable chargpd particlps of any kind with massps bplow 
100 GpV. Wp also rpquirp that any npw contributions to thp invisiblp width of thp Z 
boson bp < 2 MpV[23]; this constraint pliminatps thp possibility of certain sppcips of 
nputralinos having massps bplow M^/2. 

Following ALEPH[23] wp implpmpnt a lowpr limit of 92 GpV on first and spcond 
gpnpration squark massps, providpd that thp gluino is morp massivp than thp squarks 
and thp mass diffprpncp (Am) bptwppn thp squark and thp LSP is > 10 GpV. Wp also 
implpmpnt a similar cut (following [25]) on thp mass of sbottom quarks rpquiring that 
thpir mass bp grpatpr than 95 GpV (in addition to Am > 10 GpV, and thp mass bping 
Ipss than thp gluino mass). Thp situation for stops is slightly morp complpx|26]: wp 
dpmand that thp lightpst stop mass bp grpatpr than 97 GpV if thp stop is too light to 
dpcay into Whx\- If thp stop can dpcay to ihi); wp havp a lowpr limit of 95 GpV on its 
mass. 

Following [26], wp dpmand that right-handpd slpptons havp massps grpatpr than 100, 
95, or 90 GpV for splpctrons, smuons, and staus rpsppctivply. We only apply this limit 
whpn thp condition 0.97msiepton > ''^lsp is satisfipd. We can also apply thpsp bounds to 
Ipft-handpd slpptons, providpd that thp nputralino t— channpl diagram may bp npglpctpd 
in thp casp of splpctrons; wp assump that this is thp casp. 

Wp dpmand that chargino massps bp grpatpr than 103 GpV, providpd that thp LSP- 
chargino mass splitting is Am > 2 GpV [26]. If Am < 2 GpV, thp bound is 95 GpV. 
It should bp notpd that whpn Am is vpry small ( ^ 100 MpV), thp chargino is stablp 
for dptpctor Ipngth scalps and thp modpl will bp pxcludcd by thp stablp charged particle 
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constraints. In the case where the hghtest chargino is dominantly Wino, we can only 
apply this limit when the electron sneutrino t— channel diagram is negligible; we take 
this to be the case when the electron sneutrino is more massive than 160 GeV. 

The LEP Higgs Working Group [27]. provides five sets of constraints on the MSSM 
Higgs sector imposed by LEPII data. These are essentially limits on the Higgs-Z 
coupling times the branching fraction for decay to given final states, as a function of the 
Higgs masses. We employ SUSY-H1T[28J to analyze these. We include a theoretical 
uncertainty on the calculated mass of the lightest Higgs boson of approximately 3 
GeV[29j when applying these constraints. 

2.4- Tevatron Constraints 

We also employ constraints from the Tevatron. We obtained restrictions on the squark 
and gluino sectors arising from the null result of the DO multijet plus missing energy 
search [30]. We generalize their analysis to render it model independent, by generating 
multijet plus missing energy events for our model spectrum using PYTHIA6.4[3T] (which 
we provide with a SUSY-HIT[2H] decay table) as interfaced to PGS4 [32]. PGS4 provides 
a fast detector simulation and is used to impose the kinematic cuts used in the DO 
analysis. We weigh our results with K factors computed using PROSPINO2.0 [33]. The 
95% CL upper limit on the number of signal events, as defined by the DO analysis, is 
8.34 (for the 2.1 fb~^ data set considered) using the method of Feldman and Cousins [M]. 
Analogously, we employ constraints from the CDF search for trileptons plus missing 
energy [35], which we also generalize to the full pMSSM. We only make use of the CDF 
'3 tight lepton' analysis as it is the cleanest and easiest to implement with PGS4; we also 
use a K-factor of 1.3 for all models. Here the 95% CL upper bound on the possible SUSY 
signal in the channel we are considering is 4.65 events for the luminosity of 2.02 fb~^ 
used in the CDF analysis. 

In addition to these collider signature bounds, we also employ the experimental 
constraint [36] resulting from direct searches for the new Higgs fields in the MSSM: for the 
narrow mass range 90 < < 100, tan f] is restricted to the region tan /3 > 1.2m^ — 70. 
This range is excluded as the Tevatron would have otherwise discovered at least one of 
the heavier Higgs bosons. Also, D0[37J has obtained lower limits on the mass of heavy 
stable charged particles. We take this constraint to be > 206\Uiy^\'^ + 171|f/i/ip 
GeV at 95% CL for charginos, where the matrix entries Uiw and Uih determine the 
Wino/Higgsino content of the lightest chargino. We use this to interpolate between the 
separate Wino and Higgsino results provided by DO. 

CDF and DO also have analyses that search for light stops and sbottoms^38j which 
include a number of assumptions about the SUSY mass spectrum, sparticle decay 
channels, etc. In general they are only applicable when the sbottoms or stops are lighter 
than the top quark. These searches are difficult to implement in a model-independent 
pMSSM context. Thus we exclude models with light (m < mt) stops or sbottoms from 
our final set of models; this only affected ~ 1000 models. 
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2.5. Astrophysical Constraints 

There are two constraints from considering the LSP as a long-hved rehc. As noted 
above, we demand that the LSP be the hghtest neutrahno. We also require, following 
the 5 year WMAP measurement [39] of the relic density, that f2/i^|LSP ^ 0.121. In not 
employing a lower bound on Q^^Ilsp for our models, we acknowledge the possibility 
that even within the MSSM and the thermal relic framework, dark matter may have 
multiple components with the LSP being just one possible contributor; we thus only 
require that the LSP not have a relic density too large to be consistent with WMAP. 
However, in discussing results below, we will also discuss a subset of models for which 
0.1 < Qh'^lisp < 0.121; these represent the more standard assumption that the LSP is 
the dominant, perhaps only, component of the relic density. 

We also obtain constraints from attempts to detect dark matter directly [lO]. 
Generally, the strongest constraints come from the spin- independent WIMP-nucleon 
cross sections, hence we only implement bounds on our models from these; inspection 
of the spin-dependent WIMP-nucleon cross sections in our models confirms that this 
approach is reasonable. Both spin-independent and spin-dependent cross sections 
were calculated using micrOMEGAs2.21[9]. We implement cross section limits from 
XENONlOm], CDMSt42j, CRESST 103] and DAMAgl] data. Since these cross 
sections depend on some low energy quantities for which the uncertainties are relatively 
large {e.g., nuclear form factors), we do not exclude models with WIMP-nucleon spin- 
independent cross sections as much as 4 times larger than the experimental bounds. It 
should be noted that many of our models predict a value f2/i^|LSP which is less than that 
observed by WMAP and supernova searches. We thus scale our cross sections to take 
this into account. 

3. Results 

As noted above we randomly generated 10^ parameter space points {i.e., models) in a 
19-dimensional pMSSM parameter space using flat priors. Only ~ 68.5 • 10^ of these 
models satisfy all the constraints listed in the previous section. The properties of these 
models are described in much greater detail in [8]. Here we will discuss the attributes of 
these models which are most important astrophysically. In particular we will examine 
the mass and composition of the LSP and nLSPs, the predicted relic density, as well as 
direct and indirect dark matter detection signals from these models. 

3.1. LSP and nLSP 

Figure [1] presents a histogram of the masses of the four neutralino species in our 
models; Figure [2] displays a similar histogram for the two chargino species. The lightest 
neutralino is, of course, the LSP. The LSP mass lies between 100 and 250 GeV in over 
70% of our models. Generally models with a mostly Higgsino or Wino LSP have a 
chargino with nearly the same mass as the LSP; as sufficiently light charginos would 
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normally have been detected at LEP or the Tevatron, there are fewer models with such 
LSPs with mass < 100 GeV. 
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Figure 1. Distribution of ncutralino masses for our set of models. 



The identity of the nLSP is shown in Figure El The lightest chargino is the nLSP 
in about 78% of the models; this is due to many models having Wino or Higgsino LSPs, 
and the generally small mass splitting between a mostly Wino or Higgsino neutralino 
and the corresponding chargino. The second lightest neutralino is the nLSP ~ 6% of 
the time. These will generally be models with a dominantly Higgsino LSP. Note also 
that while neutralinos or charginos are the nLSP in the vast majority of cases, there 
are 10 other sparticles each of which is the nLSP in > 1% of our models. Scenarios in 
which these sparticles are the nLSP may lead to interesting signatures at the LHC|45j. 

Figure m displays the LSP mass value as a function of the LSP-nLSP mass splitting. 
Am, our models for each identity of the LSP. It is interesting that these models have a 
smaller Am than is often considered; 80% of our models have Am < 10 GeV, 27% have 
Am < 1 GeV, and 3% have Am < 10 MeV. As one can see from Figure HI this occurs 
largely, but not exclusively, in models with a chargino nLSP. This is again due to the 
many models where the LSP is nearly pure Wino or Higgsino. 

There are a number of interesting features in this figure. The mostly empty square 
region which appears on the lower left-hand side of Figure H] is due to the fact that 
models with chargino nLSPs in this mass and Am range have been excluded by the 
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Figure 2. Distribution of chargino masses for our set of models. 



Tevatron stable chargino search. Non-chargino nLSPs are not ehminated by this search 
{e.g., the production cross section for sleptons in this range is too small to be excluded 
by the Tevatron search). It is perhaps worth noting that a stable heavy charged particle 
search at the LHC, corresponding to those done at the Tevatron, would be able to 
exclude or discover the models with heavier chargino nLSPs and small values of Am 
(corresponding to ~ 12% of our model set). 

Another interesting feature in this figure is the bulge for 0.1 GeV < Am < 2 GeV 
and rriLSP ~ 100 GeV. This region exists because these values of Am are large enough 
that at LEP or the Tevatron, the produced chargino would decay in the detector, but 
the resulting charged tracks would be too soft to be observed. The existence of such 
a region shows the difficulty of making model independent statements about sparticle 
masses or other SUSY observables. 

We have seen that within our model set the nLSP can be almost any SUSY particle 
and the corresponding Am can be small for these cases. Thus specific models in our 
set describe qualitatively most of the conventional long-lived sparticle scenarios. Long- 
lived stops or staus (as in GMSB[i2j), gluinos (as in Split SUSY [16]) as well as charginos 
(as in AMSB[3]) all occur in our sample. We also have long-lived neutralinos, as does 
GMSB, however these are the X2 ou'^ case. In addition to models which, to some 
extent, correspond to these well-studied scenarios, we also have models with long-lived 
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Figure 3. Number of models in which the nLSP is the given sparticle. 

selections, sneutrinos and sbottoms. 

Figures O O and [7] display the gauge eigenstate content of the LSPs in our model 
set. We note that most LSPs are relatively pure eigenstates, with models where the 
LSP is Higgsino or mostly Higgsino being by far the most common. About one quarter 
of our models have Wino or mostly Wino LSPs, while just over one-sixth have Bino or 
mostly Bino LSPs. Within mSUGRA, the LSP is, in general, nearly purely Bino; this 
suggests that most of our models are substantially different from mSUGRA. A more 
precise breakdown of the content of LSPs in the model set is presented in Table 
We note that one would expect the LSP be a pure eigenstate fairly often in a random 
scan of Lagrangian parameters, since if the differences between Mi, M2, and are large 
compared to M^, then the eigenstates of the mixing matrix will be essentially pure 
gaugino and Higgsino states 0. 

3.2. Relic Density 

We did not demand that the LSP, in any given model, account for all of the dark matter, 
rather we required only that the LSP relic density not be too large to be consistent 
with WMAP. More specifically, we employed ^/i^Ilsp < 0.121. Figure [H] shows the 
distribution of ^/i^Ilsp values predicted by our model set. Note that this distribution 
is peaked at small values of f2/i^|LSP- In particular, the mean value for this quantity in 
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Figure 4. Mass splitting between nLSP and LSP versus LSP mass. The identity of 
the nLSP is shown as well. (The LSP is always the lightest neutralino in our set of 
models). 



LSP Type 


Definition 


Fraction 
of Models 


Bino 


|Zn|2 > 0.95 


0.14 


Mostly Bino 


0.8 < |Zii|2 < 0.95 


0.03 


Wino 


|^i2p > 0.95 


0.14 


Mostly Wino 


0.8 < \Zi2? < 0.95 


0.09 


Higgsino 


1^131' + I^Mp > 0.95 


0.32 


Mostly Higgsino 


0.8 < |Zi3|2 + I^Mp < 0.95 


0.12 


All other models 




0.15 



Table 2. The fractions of our model set for which the LSP is of each of the given 
types. These types are defined here by the modulus squared of elements of neutralino 
mixing matrix in the SLHA convention. See |5j for details. 
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Figure 5. The distribution of LSP gaugino eigenstate types as a function of the LSP 
mass. Note that each LSP corresponds to three points on this figure, one each for its 
Bino, Wino, and Higgsino fraction. 



our models is ~ 0.012, which is about ten times less than the central value determined 
from the aforementioned WMAP and supernova data[39j. We note that the range of 
possible values of O/i^Ilsp is found to be much larger than those obtained by analyses 
of specific SUSY breaking scenariosPTj. 

We display the predictions for f2/i^|Lgp versus the LSP mass in Figure [9] and 
versus the nLSP - LSP mass splitting in Figure dUJ Figure [HI makes it clear that 
f2/i^|LSP generally increases with the LSP mass, but a large range of values for the 
relic density are possible at any given LSP mass. The empty region in Figure [9] where 

LSP ~ 0.001 — 0.1 and mLSP ~ 50 — 100 is due to the fact that, in general, LSPs 
which are mostly Higgsino or Wino give lower values of VLh?\i^^^, and, as noted above, 
there are fewer Higgsino or Wino LSPs in this mass range. Figure [ID] shows that small 
mass differences can lead to large dark matter annihilation rates. 

3. 3. Direct Detection of Dark Matter 

As noted above, we calculate the spin- dependent and spin-independent WIMP-nucleon 
cross sections using micrOMEGAs 2.21 [9]. These data give the possible signatures in 
our model set for experiments that search for WIMPs directly. As these experiments 
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nLSP-LSP 

Figure 6. The distribution of LSP gaugino eigenstate types as a function of the LSP- 
nLSP mass difference. Note that as in Figure [51 every LSP corresponds to three points 
on this figure, one each for its Bino, Wino, and Higgsino fractions. 

measure the product of WIMP-nucleon cross sections with the local relic density, the 
cross section data presented in the figures below are scaled by ,^ = ^/i^Ilsp/^^^Iwiviap- 
To date, these experiments generally provide a more significant bound on the spin- 
independent cross section, and hence we will focus on those. 

Figure [TT] presents the distribution for the scaled WIMP-proton spin-independent 
cross section versus relic density for our model sample. As one would expect, larger 
values of the cross section are generally found at larger values of Q/i^Ilsp. However, 
even for rehc densities close to the WMAP value, ^(Jp^si is seen to vary by almost eight 
orders of magnitude. These ranges for ^(Tp^si are much larger than those from mSUGRA 
as calculated, e.g., in [18] . 

Figure [12] shows the scaled WIMP-proton spin-dependent and spin-independent 
cross sections as a function of the LSP mass. The constraints from XENON10j41j and 
CDMS[l2] are also displayed. As noted above, to take the uncertainties in the theoretical 
calculations of the WIMP-nucleon cross section into account, we allowed for a factor 
of 4 uncertainty in the calculation of the WIMP-nucleon cross section. Table [3] gives 
the fraction of models that would be excluded if the combined CDMS/XENONIO cross 
section limit were improved by an overall scaling factor. Note that our inclusion of the 
theoretical uncertainties does not significantly modify the size of our model sample. 
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Figure 7. Wino/Higgsino/Bino content of the LSP in the case of flat priors. Note that, 
as elsewhere in the paper, |2'i2p, and \Zi^\'^ + |Zi4p, where Zij is the neutrahno 

mixing matrix in the SLHA convention [S], give the Bino, Wino, and Higgsino fractions 
respectively. 
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Figure 8. Distribution of fi/i^|LSP for our models. 



LSP Mass Versus Relic Density 
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Figure 9. Distribution of ^h^\'Lsv as a function of the LSP mass. 
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LSP - nLSP Mass Difference Versus Relic Density 
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Figure 10. Distribution of JI/i^Ilsp as a function of the LSP-nLSP mass splitting. 



Improvement in S.I. 
Cross Section Limit 


Fraction of Models 
Excluded 


4 


0.032 


10 


0.071 


40 


0.19 


100 


0.31 


400 


0.52 


1000 


0.65 


4000 


0.81 



Table 3. The fraction of our model set which would be excluded for the specified 
improvement in the direct detection bound on the spin-independent WIMP-nucleon 
cross section. 

We find that the range of values obtained for these cross sections covers the entire 
region in cross section/ LSP space that is anticipated from different types of Beyond the 
Standard Model theories in the above reference. This possibly suggests that we cannot 
use direct detection experiments to distinguish between e.g. SUSY versus Little Higgs 
versus Universal Extra Dimensions dark matter candidates in the absence of other data. 

In Figure [131 we compare the WIMP-proton and WIMP-neutron cross sections in 
the spin-dependent and spin-independent cases. The spin-independent cross sections 
are seen to be fairly isospin independent; this is not the case, however, for the spin- 
dependent cross sections. 
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WIMP-Prolon ag. Versus Relic Density 




1e-18 I • • — -I ■ ^ • ^ • • — -I ■ ^ ■ ^ 

1e-06 1e-05 0.0001 0.001 0.01 0.1 1 

Figure 11. Distribution of scaled WIMP-proton spin-independent cross section versus 
the LSP contribution to relic density for our models. 

3.4- Indirect Detection of Dark Matter 

The PAMELA collaboration has recently claimed an excess in the ratio of cosmic ray 
positrons to electrons observed at energies > 10 GeV[49j. Here we employ DarkSUSY 
5.0.4[50] to calculate this ratio for our model sample and compare these results with the 
PAMELA data. 

In general, for a thermal relic dark matter candidate to reproduce the PAMELA 
data, its signal rate must be multiplied by a boost factor [51]. In nature, such a boost 
factor could result from, e.g., a local overdensity. The boost factor in that case would 
be the square of the ratio between the density of dark matter in the region from which 
one is sensitive to cosmic ray positrons and electrons to the universe as a whole. 

We have investigated four of the propagation models available as default choices 
in DarkSUSY: the model of Baltz and Edsjo(BE)^2], that of Kamionkowski and 
Turner (KT) [53], that of Moskalenko and Strong(MS)[53], as well as GALPROP[55]. 
In the figures that follow we show the results of calculations using the MS propagation 
model, which is based on early GALPROP Green's functions and whose results typically 
lie between those computed otherwise. However, we note that the extent to which 
the positron/electron flux ratio predicted by our models matches the PAMELA data 
can be highly sensitive to the choice of propagation model parameters and assumed 
astrophysical backgrounds. We will explore this further in future work|56j. The halo 
model employed here is the Navarro- Frenck- White profile [57]. 

The differential positron flux as a function of energy for a random sample of 500 
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Figure 12. Distributions of scaled WIMP-proton spin-dependent cross section 
and spin-independent cross sections versus LSP mass in our models. In the spin- 
independent panel, the CDMS and XenonlO bounds, which provide the strongest limits 
for the range in LSP mass relevant for our models, are shown. 
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Figure 13. Here we compare WIMP-neutron and WIMP-proton cross sections. The 
spin-dependent cross sections are shown in the top panel; the spin-independent cross 
sections in the bottom paneL 



Dark Matter m the MSSM 



18 




Figure 14. Expected, scaled, contribution to the absolute flux of positrons 
(unboosted) from neutralino annihilation in the halo for 500 randomly selected models 
employing MS propagation. Also shown are data from HEAT^GIJ, AMS01[62]. and 
CAPRICE94[63] as well as the DarkSUSY default secondary positron background (a 
parameterization of model 08-005 in|60j. 
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Figure 15. The distribution of per degree of freedom versus the choice of boost 
factor that minimized this quantity for 500 randomly selected pMSSM models in our 
model set for the four propagation models discussed in the text. 
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Figure 16. Positron/ electron flux ratio versus energy taking models from a set of 500 
random pMSSM models for which the per degree of freedom with the x^-maximizing 
boost was less than 10.0 for three of the four propagation models studied (curves shown 
use MS propagation). Data shown are from HEAT[64], CAPRICE94[63], PAMELA[49] 
and AMS01[62]. 

models from our set are shown in Figure [TH Here we assume a boost factor of 1; the 
normahzation of the curves takes into account the fact that for many of these models 

^^^IlSP < ^WMAP- 

We next determine how well the predicted positron fluxes for these models agree 
with the PAMELA data, allowing for the possibility of a boost factor. To do this, we find 
the value for the boost factor (with the restriction that it be < 2000) which minimizes 
the foi' the fit of each model's prediction to the PAMELA data. (Note that many 
of the models require an even larger boost to obtain a good fit and are thus not shown 
in Figure [T5l) . In calculating the x^? we consider only the seven highest energy bins, 
as at lower energies solar modulation is expected to play a major role[l9]. Figure fT5l 
shows the ^ind corresponding boost factor for these 500 random models. Note that 
there are four data points for each model in this figure. We then display the positron 
to electron flux ratio, for the models with a low value of x^ employing MS propagation, 
as a function of energy in Figure [161 and note the reasonable agreement with the data 
for some models. 

Since the flux from WIMP annihilation scales as (f2/i^|LSp/^^^|wMAp)^, we might 
expect to improve the match to the PAMELA data using models from our sample for 
which ^/i^Ilsp ~ ^^^IwMAP- To test this, we examine the predicted positron flux for 
500 random models with i7/i^|LSP > 0.100 employing MS propagation; these fluxes are 
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Figure 17. Same as Fig. 14 but now for a 500 model set satisfying ^h?\wMAP > 
f^/i^ksp > 0.10. 




Figure 18. Same as Fig. 15 but now for a 500 model set satisfying il/i^lwMAP > 
nh\sv > 0.10. 
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Figure 19. Positron/ electron flux ratio versus energy taking models from a set of 500 



pMSSM models satisfying Q.h^ 



WMAP 



> rj/i I LSP > 0.10 and for which x maximizing 



boost was less than 5.0 for three of the four propagation models studied (curves shown 
use MS propagation). Data shown are from HEAT[64], CAPRICE94[63], PAMELA[49] 
and AMSOlIMl. 



shown in Figure [T7] with no boost factor. We then again find the boost factor that 
minimizes the of the positron to electron fiux ratios with respect to the seven highest 
energy PAMELA bins; these are shown in Figure [HI Here, we note that there are many 
more models for which the x^-minimizing value for the boost factor is < 2000 and there 
are many more points for which the value is low. The positron to electron fiux ratios 
for these models, including the boost factor, are shown in Figure [191 

It appears that some of our models do a reasonably good job of fitting the PAMELA 
positron data, especially in the case where ^/i^Ilsp lies fairly close to the WMAP value. 
For most models, describing the PAMELA data requires large boost factors, however 
this is also a fairly generic feature of attempts to explain PAMELA and ATIC data in 
terms of WIMP annihilation [5 Ij. There are however, many models which give relatively 
low per degree of freedom in the fit to the data with relatively small boost factors. 
We will study this further in future work [56]. A study of the corresponding predictions 
for the the cosmic ray anti-proton fiux is also underway. 

4. Conclusions 

We have generated a large set of points in parameter space (which we call "models" ) for 
the 19-parameter CP-conserving pMSSM, where MFV has been assumed. We subjected 
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these models to numerous experimental and theoretical constraints to obtain a set of 
~ 68 K models which are consistent with existing data. We attempted to be somewhat 
conservative in our implementation of these constraints; in particular we only demanded 
that the relic density of the LSP not be greater than the measured value of VLh"^ for non- 
baryonic dark matter, rather than assuming that the LSP must account for the entire 
observed relic density. 

Examining the properties of the neutralinos in these models, we find that many are 
relatively pure gauge eigenstates with Higgsinos being the most common, followed by 
Winos. The relative prevalence of Higgsino and Wino LSPs leads many of our models 
to have a chargino as nLSP, often with a relatively small mass splitting between this 
nLSP and the LSP; this has important consequences in both collider and astroparticle 
phenomenology. 

We find that, in general, the LSP in our models provides a relatively small (~ 4%) 
contribution to the dark matter, however there is a long tail to this distribution and 
a substantial number of models for which the LSP makes up all or most of the dark 
matter. Typically these neutralinos are mostly Binos. 

Examining the signatures of our models in direct and indirect dark matter detection 
experiments, we find a wide range of signatures for both cases. In particular, we find, 
not unexpectedly a much larger range of WIMP-nucleon cross sections than is found 
in any particular model of SUSY-breaking as can be seen by comparing directly with 
the work in Ref . [48] . In fact, as these cross sections also enter the regions of parameter 
space suggested by non-SUSY models, it appears that the discovery of WIMPs in direct 
detection experiments might not be sufficient to determine the correct model of the 
underlying physics. As a first look at the signatures of these models in indirect detection 
experiments, we examined whether our models could explain the PAMELA excess in 
the positron to electron ratio at high energies. We find that there are models which fit 
the PAMELA data rather well where the LSP is mostly Bino, and some of these have 
significantly smaller boost factors than generally assumed for a thermal relic. 

The study of the pMSSM presents exciting new possibilities for SUSY 
phenomenology. The next few years will hopefully see important discoveries both in 
colliders and in satellite or ground-based astrophysical experiments. It is important 
that we follow the data and not our existing prejudices; hopefully this sort of relatively 
model-independent approach to collider and astrophysical phenomenology can be useful 
in this regard. 
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